Decoherence-free quantum information in the presence of dynamical evolution 
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We analyze decoherence-free (DF) quantum information in the presence of an arbitrary non- 
nearest-neighbor bath-induced system Hamiltonian using a Markovian master equation. We show 
that the most appropriate encoding for N qubits is probably contained within the ~ |JV excitation 
subspace. We give a timescale over which one would expect to apply other methods to correct 
for the system Hamiltonian. In order to remain applicable to experiment, we then focus on small 
systems, and present examples of DF quantum information for three and four qubits. We give an 
encoding for four qubits that, while quantum information remains in the two-excitation subspace, 
protects against an arbitrary bath-induced system Hamiltonian. Although our results are general 
to any system of qubits that satisfies our assumptions, throughout the paper we use dipole-coupled 
qubits as an example physical system. 



I. INTRODUCTION 

Storing quantum information for long periods is dif- 
ficult: excited quantum states decohere. One method 
of counteracting the effects of decoherence is to encode 
logical information into decoherence-free subspaces and 
subsystems (DFSs) 0, H, i, 1, i, i, 0, @. These are 
groups of states that have robust symmetry properties, 
and that in ideal cases serve as perfect quantum mem- 
ory. See Ref. 0] for a comprehensive review of DFS the- 
ory. As well as theoretical developments, there has been 
much success experimentally. For example, DFSs have 
been prep ared in optical systems, in NMR, and in ion- 
traps [ljj [HI, [lH, [l]| , and there has been a proposal for 
decoherence-free (DF) quantum-information pro cessing 
in nitrogen- vacancy (NV) centers in diamond [141 ]. 

Here, we focus on the semigroup formulation of DFSs 
for a number of reasons. First, the unitary and nonuni- 
tary evolution of the system are naturally separated. Sec- 
ond, when the DF condition is satisfied, there is only 
one nonunitary Lindblad operator 0, Q. Finally, the 
scalability properties of various encodings can be quanti- 
fied @, ■ The separation of unitary and nonunitary evo- 
lution enables the detrimental effect of any environment- 
induced unitary evolution to be isolated. We refer to 
'environment-induced', or equivalently 'bath-induced', 
evolution because both the unitary and nonunitary parts 
of the master equation result directly from the system- 
environment interaction Hamiltonian. For a single atom, 
the unitary part gives rise to the Lamb-shift and the 
nonunitary part gives rise to spontaneous decay. For 
many qubits, decoherence can result indirectly from the 
action of Lamb-shift type terms because these can cause 
quantum information to leak into states that decay. A 
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full analysis of the effect of unitary evolution on DF quan- 
tum information in Markovian systems is the purpose of 
this paper. Throughout the paper, when we refer to a 
system Hamiltonian causing leakage of DF quantum in- 
formation, we are referring explicitly to the action of a 
Lamb-shift type Hamiltonian on DF states. 

It is well-known that the condition for DF dynamics 
derived in Ref. Q does not rule out the possibility of 
an environment-induced unitary operator evolving infor- 
mation encoded in a DFS into non-DFS states. The re- 
sults in Ref. Q were extended in Ref. [l5| , where the au- 
thors gave conditions for DF quantum information that 
accounted for a system Hamiltonian. Here, we focus on 
quantum information that is encoded within states that 
satisfy the DF condition 0, but that could evolve into 
non-DF states due to a system Hamiltonian. We em- 
phasize that for Markovian systems, the property of the 
bath operators is such that generally, the presence of 
a non-nearest-neighbor system Hamiltonian is unavoid- 
able. Physical justification for the regime studied in this 
paper is given in Ref. [l6| . Here, the authors showed that 
for closely-spaced dipole-coupled qubits, the nonunitary 
evolution can be approximated as being the same as that 
for co-located qubits (and so satisfy the DF requirements 
derived in Ref. Q), but the unitary evolution depends on 
both the position and orientation of the dipoles. Our re- 
sults are applicable to any physical system that can be 
appropriately described by a Markovian master equation 
and, in light of recent experimental progress, are directly 
applicable to NV centers in diamond [17l ]. 

The paper is structured as follows. In Sec. [H] we state 
the conditions on the bath operators that lead to a bath- 
induced unitary evolution. In Sec. IIII1 we give a con- 
dition for subspaces which are immune to nonunitary 
evolution. This condition is weaker than that derived 
in Ref. but stronger than that derived in Ref. 0]. 
Then, in Sec. IIVI we analyze the scalability — which we 
define as the encoding efficiency multiplied by the pro- 
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portion of the Hilbcrt space that is DF — for DF encoding 
in the presence of a system Hamiltonian. An interesting 
consequence of our results is that, when the effect of a 
system Hamiltonian is included, the most suitable sub- 
space for quantum information storage in N qubits is 
probably the subspace with ~ |iV excitations. Then, we 
give a timescale over which other methods will have to 
be applied to counter the effect of the system Hamilto- 
nian. In Sec. |Vl we concentrate on three and four qubit 
systems. We give conditions on the system Hamiltonian 
for DF quantum information, and we give an explicit en- 
coding for four qubits that, while quantum information 
remains in the two-excitation subspace, protects against 
errors induced by an arbitrary system Hamiltonian. We 
hope that this encoding is directly relevant to experimen- 
tal quantum information processing. 

II. MASTER EQUATION 

For clarity and to establish our notation, we briefly 
summarize a derivation of the Lindblad master equation. 
A system A coupled to a bath B can be described by 
the Hamiltonian H = <£> Ig + Ia <8 Hb + H/, where 
H_4(Hg) the system (bath) Hamiltonian acts on the sys- 
tem (bath) Hilbert space, I^(Ib) is the identity operator 
on the system (bath) Hilbert space, and H/ is the inter- 
action Hamiltonian that contains all non-trivial couplings 
between the system and the bath and is written as 

H/ = ^S Q ®B a , (1) 

OL 

for S a (B a ) the system (bath) operators. Treating the 
interaction Hamiltonian as a perturbation, the equation 
of motion for the density matrix x of the system and bath 
in the interaction picture is (H = 1) 

X (t) = -i[tt I (t), X (t)] 1 (2) 

which gives 

p(t) = -[ dsTz B [H I (t),[li I (s),x(s)]}, (3) 
Jo 

where Tr^ denotes the trace over the bath and p is the 
(reduced) system density matrix. So, for the rest of this 
paper, we focus on the Hilbcrt space of the system. Mak- 
ing the Born-Markov approximations, Eq. ^ becomes 

/>oo 

p(t) = - dsTr B [H J (t),[H J (s),p(t)®i? ]] (4) 

for Rq the initial density operator of the (stationary) 
bath. Introducing the correlation function 

C a ^s) = Ti B [Bi(s)B (O)Ro}, (5) 

which we note has real and imaginary parts because in 
general B^(s) and B^(0) do not commute, the mas- 



ter equation in the rotating-wave approximation is writ- 
ten (l8j 

p = L[p}=-i[n s ,p]+~L B [p], 
l d[p] = ^2a a pL SatSfj [p], 

a, /3 

Ls B1 s>] = [S /Jl pSt] + [S /3 p,St] 1 (6) 
for the system Hamiltonian 

H S = J^St S,,, ( 7 ) 

and where a a p = T a/3 + T* Q and b al3 = ± [T a p - r* J for 
r Q/3 = / dse^C^s). (8) 

The correlation function C a p(s) contains all the in- 
formation about the physical system, and satisfies the 
Kramers-Kronig relations. For a single-atom, Hs de- 
scribes the Lamb-shift, and Ld[p] describes the sponta- 
neous emission. For two or more dipole-coupled qubits, 
the off-diagonal terms in the Hermitian matrix (a a p) de- 
scribe the rate of spontaneous emission between separate 
physical qubits, and the off-diagonal terms in the ma- 
trix (b a p) describe the coherent dipolc-dipolc interaction 
between separate qubits. Generally, the terms in (b a p) 
diverge and require renormalization [l6j. 

A. Decoherence-free subspaces 

A subspace p of the system Hilbert space p is a DFS 
if it satisfies Ld[/5] =0 [3]. The condition that leads to 
Ld[p] = being satisfied is a a p = a [Hj]. For dipole- 
coupled qubits, this implies that the spontaneous emis- 
sion rate between separated qubits is the same as that for 
individual qubits — all qubits experience the same nonuni- 
tary couplings. Here, we focus on the effect of the sys- 
tem Hamiltonian on p. Note that when we refer to the 
system Hamiltonian, we are not referring to H^, but 
Hs which results from H/. We consider the situation 
in which the Hamiltonian Hg evolves information en- 
coded in a DFS into non-DFS states @, [20j]. We are in- 
terested in the physically applicable perturbative regime 
(o Q /3) = aX + eA and (bap) = 6X + eB for X a matrix 
of size (a a p) with all entries equal to one, A and B ma- 
trices also of size (a a p), but with arbitrary entries, and 
expansion parameter e«l. 

Specifically, we concentrate on the regime (a a p) = aX 
and (b a p) = 6X + eB. We do this for two reasons. First, 
it has been shown that DFSs are stable to first order even 
in the presence of a nonunitary symmetry breaking per- 
turbation [2(|. All else being equal, small changes in the 
form of Ld[/c] do not prevent infinite-lifetime quantum 
information storage. Second, the strict requirement that 
all qubits experience the same environment is unlikely 
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to ever be met in practice. If one perturbs the system- 
environment couplings, the perturbation applies to both 
the nonunitary and unitary parts of the master equation. 

We illustrate the regime of interest to this paper using 
dipole-coupled qubits. For these qubits the appropriate 
expansion parameter is physical separation r. So, we 
concentrate on collections of qubits that satisfy 



a a /3 = a + C(r 2 ), 

b a p = b + b' ap (r)+0(r 2 ), 



(9) 
(10) 



where the system of qubits is such that any 0(r) con- 
tribution to a a p can be neglected. Eqs. ((9]) and (fT0|) 
are satisfied by closely-spaced dipole-coupled qubits — 
explicit forms for a a p and b a p are given in Ref . [1 61 ] . Note 
that b' a p{r) can be many orders of magnitude larger than 
a. 

We emphasize that the results presented here are ap- 
plicable to any Markovian system for which A = and 
B ^ 0. This is the most general case for Markovian sys- 
tems that satisfy the DF condition, and is more likely to 
be realised in the laboratory than A = and B = 0. 
Note that if A = and B = 0, then T a p is strictly 
qubit independent, and does not evolve information 
encoded in p into non-DF states 0, @ . 

For non-Markovian systems, leakage resulting from 
Lamb-shift type terms in the master equation can be 
accounted for using dynamical decoupling methods, or 
'bang-bang' pulses [U, [z| [H, |24[- Here, we focus on 
passive error correction that does not require fast and 
strong pulses and so is more amenable to experimen- 
tal implementations. Also, bang-bang pulses are inher- 
ently non-Markovian, and applying similar techniques in 
a Markovian environment is difficult. 



III. COMPLETELY-DECOHERENCE-FREE 
SUBSPACES 



The condition a a p 
be written 



a means that the dissipator can 



L D [p] = A(2JpJ f - J f Jp - pjtj), 



(11) 



where A is the only nonzero eigenvalue of a a p and J = 
S Q for N the number of qubits. The jump operator 
can be written 



N 



N 



E*<- = E" 



<7_ g)J ■ ■ ■ 
i th qubit 



(12) 



where / is the 2x2 identity operator. So, &i a , for a = 
+,—,z, acts on the i th qubit, and satisfies [<Ji Z ,6'j±} = 
±25ij&i±, and = Sij&iz- The operators J, J', 

and J 2 act on the system Hilbert space TL = ® N C 2 . From 
the representation theory of su(2), H. can be decomposed 
into irreducible components [25[ — see Fig [3] for a three 
qubit example. Each irreducible representation (irrep) V 



is generated by a unique lowest weight vector v satisfying 
J • u = J_ • v~ = 0. The vector v~ is an eigenvector of 



N 



E 5 * 



(13) 



to |e). 

—a subspace W of 



The vector J fc • v, for k applications of J on v is an eigen- 
vector of J z with eigenvalue — (N — 2k). The decoher- 
ence operator J causes the system to decay to its lowest 
weight. The DF condition Ld[V] = for V a subspace 
implies that K- V = for all jump operators K Q- Here, 
we only have one jump operator J, so the subspace V sat- 
isfies J • V = 0, and consists of combinations of lowest 
weight vectors. For example, the states |b) and |c) in 
Fig [3] satisfy J|b) = and J |c) = 0. They are states 
that do not decay to |a) and so are DF. The operator J' 
acts like Jt |b) = |e) and Jt |c) = |f), but does not cause 
transitions from |b) to |f), or from 
The definition for DF dynamics 
states is DF if any p(0) G W evolves into a state p(t) 
such that the evolution map p(0) — ► p(i) is unitary V t — 
is relevant here. We are interested in the evolution of 
quantum information under the action of both Ld and 
Hg. Ensuring Ld[p] = for some subspace is not suf- 
ficient to guarantee DF dynamics for all time. We call 
a subspace W completely-decoherence-free (CDF) if it 
satisfies the following conditions 

1. L D [W] = and 

2. p(t)eW Vt. 

These conditions ensure that the evolution is unitary. We 
note that the result at the replica symmetric point de- 
rived in Refs. [J, Q ensures CDF dynamics, and that 
the conditions stated here still permit the transfer of 
encoded information between states within W . For the 
purposes of this paper, when we refer to transfer of en- 
coded information, we are referring to the effect of the 
off-diagonal terms of Us in the Clebsch-Gordan basis 
on quantum information. A basic characterization of 
completely-decoherence-free subspaces is given below. 

Proposition 1 Let V be the subspace of lowest weight 
vectors. A necessary and sufficient condition that V con- 
tain a CDF subspace W is Hg • W C W . In particular, 
Hg can be diagonalized in W . 

Proof: Define p'{t) = e lHst p(t)e~ lHst for some Hamil- 
tonian H5. The equation satisfied by p' is 



<V t / r ,1 

37 = LdH 



(14) 



where L' u [3] = L D [e lHst J e - lHst ]. So, in this pic- 
ture p'(0) is DF iff Ld[p'(0)] = 0. The generic 
DFS are spanned by vectors \x) such that J' \x) = 
e iH s tj e -jH s t j,^ _ q \y be the subspace consist- 

ing of all such vectors. This must be satisfied for all t, so 
H5 ■ W C W. Conversely, if this is satisfied then 



J ( 



-m s t 



0. 



(15) 
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□ 

This condition is weaker than that derived in Ref. [l5[, 
but stronger than that derived in Ref. Q- Throughout 
this paper, it is implicitly assumed that if Eq. (|15[) is 
satisfied, then other techniques, such as quantum error 
correction, are applied to account for transfer of encoded 
information between states within a CDFS. 

We have observed that DFSs are lowest weight vectors 
in the decomposition of C 2 into irrcps of su(2). The 
weight space W(k) of weight k in H is the eigenspacc 
of J z with eigenvalue k. Since [J z ,Hs] = weight 
spaces W(k) are left invariant by Hs- Hence, it is suffi- 
cient to consider subspaces with a fixed weight, that is, 
all states in the subspace have the same number of ex- 
cited (physical) qubits. Note that these nontrivial sub- 
spaces span across several irrcps, so it is necessary to 
combine irrcps if one is to satisfy Eq. (|15p . From the 
representation theory of su(2) [25|, it follows that the 
weight spaces have weights —N,~(N — 2), . . . , N. For 
k < N/2 the dimension of space W(k) for k excita- 
tions with weig ht —(N— 2k) is ( N k ). The condition 
J • V(k) = 0, V(k) C W(k) allows us to write the di- 
mensions of V(k) as 



dim[V{k)} 



N 
k - 1 



Nl(N-2k + l) 



k\(N 



1)1 



(16) 



IV. SCALABILITY PROPERTIES AND 
LEAKAGE TIMESCALE 

A sufficient condition for CDF dynamics derived in 
Ref. [ij for N qubits that are prepared in a DFS is as 
follows. If the unitary coupling satisfies b a p = b, then 
the Hamiltonian can be written 



Hs = &$Z S a S /' = fcjtj ' 



(17) 



a, 13 



Thus, Hg is a product of the total operators J and J' , 
and the irreps are left invariant by H5. An example 
of this is TV co-located qubits that are described by the 
master equation derived in Ref. [l6l |. In this instance, the 
multi-qubit level shift is the same for all qubits. Note 
that if b a p = b is satisfied, then DFSs are stable to a 
symmetry breaking perturbation [20| . 

For iV qubits, one can estimate the effect of Hs on en- 
coded quantum information by examining the proportion 
of the Hilbert space that consists of DF states relative to 
the proportion that consists of non-DF states. We con- 
sider vectors of the same weight, or equivalently, the same 
level of excitations. For V(k) the DF subspace in W{k) 
(with weight — (N— 2k), 2k < N) the encoding efficiency 
is defined as the number of logical qubits per number of 
physical qubits. So, using Eq. ()16|) we define 



1 



1 



dnF = — log 2 dim[V(A;)] = — log. 



N 



N\(N- 2k + 1) 
fc!(iV-fc + l)! 



(18) 




f/ D F X PDF 




(b) 



FIG. 1: (a) Plot of doF (dotted line) and pdf (solid line) 
and (b) dr>F X pdf for large N for duF and pdf as denned 
in Eqs. (|19|l and (|21[) respectively. In (b), duF X pdf is a 
maximum for r ~ §. 



The encoding efficiency c?df measures how many DF 
qubits can be encoded into a Hilbert space Tt = <g) N C 2 , 
and is unity for scalable encoding. Writing k = rN where 
r < 1/2, and taking the limit N — > 00 gives 



N- 



-rlog 2 r- (1-r) log 2 (l -r), 



(19) 



where N is the number of physical qubits, and r is in- 
dependent of N. For N — ► 00, c?df is a maximum for 
r = 1/2. This is the canonical strong-collective DFS [26[. 

Encoding efficiency alone is not sufficient as a measure 
of scalability when the Hamiltonian Hs can cause leak- 
age of quantum information from DF to non-DF states. 
So, as a measure of the likelihood that an arbitrary Hs 
causes quantum information to transfer between states 
within a particular excitation, we define 



PDF 



dim[^(jfe)] 
dim[VF(fc)] 



1 



k — 1 — N 



(20) 



which is simply the fraction of a particular weight space 
that satisfies J • V(k) = 0. So, pdf measures the propor- 
tion of DF states relative to non-DF states for a partic- 
ular excitation. For N — > 00, pdf is 



PDF 



1 - 2r 



(21) 



The smaller puf, the more likely that Hs will cause 
quantum information to evolve into non-DF states in that 
particular weight space. We show <1df and pdf for large 
TV in Fig. l(a)| and the product d^F X Pdf in Fig. 1(b) 
For an arbitrary Hs that causes quantum information 
to transfer between states within some weight space, the 
subspace which maximises djjp x pdf is one for which 
k ~ |iV. Of course, this assumes that one can account for 
the unitary evolution caused by Hs within the CDFSs. 
For particular forms of Hs, it might be more appropriate 
to encode in other weight subspaces, but if scalability in 
this context is important, then care must be taken to 
ensure djjp x puf is large. 

As an example, we focus on the case of strong-collective 
decoherence [26|. The DF subspace in this instance has 
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FIG. 2: Ratio p DP , Jtot for (i) J = 1, . . 
(iii) J = 1 for iV qubits, where PDF,j to 



AT/2, (ii) J =1,2, and 
is defined in Eq. ([23} . 



dimension 0, 0, [9[ 

dim[DFS(A0] 



AH 



(JV/2 + 1)!(A72)! : 



(22) 



for the collective basis \ J,mj), where |0)(|1)) represents 
a \j = \,fTij = — h)(\j = nij = 5)) state. Note that we 
are not referring here to the physical angular momentum 
of a particle, as in Ref. [26[ we are simply using the no- 
tation for convenience. For large N, the encoding effi- 
ciency g?df is asymptotically unity. However, for J = 
and r = 1/2, pdf N -^3° 0. So, the proportion of DF 
states in the A^/2 subspace is asymptotically zero, imply- 
ing that an arbitrary Hg that causes information to leak 
from DF-states to non-DF states negates the encoding- 
efficiency of strong-collective DF subspaces. We consider 
three cases: (i) Hs causes DF information to leak into 
all other states in the N/2 subspace, (ii) Hs causes DF 
information to leak from J = to J = 1 and J = 2 
states, and (iii) Hs causes DF information to leak from 
J = to J = 1 states. We define 



dim[DFS(A0] 
dim[non-DFS(A0] : 



(23) 



for 



dim[non-DFS(A0] = ^ 



(2J + 1)JV! 



-J (N/2 + J + l)\(N/2 — J)\ 

(24) 

Note that the quantity pr>F,j ta t equals unity for equal 
amounts of DF and non-DF states, whereas for the 
same instance jjdf equals 1/2. Fig. [2] shows pDF,j tot for 
Jtot = 1, Jtot = 2, and J tot = A72 for 500 qubits. Al- 
lowing quantum information to transfer to just one other 
subspace reduces the encoding efficiency. So, understand- 
ing the effect of H5 in particular physical realizations is 
important for scalability. 

The effect of variations away from b a p = b can be 
quantified as follows. The fidelity F(t) = Tr[pu (t) p(t)] 
for pu(t) the unwanted unitary evolution and p(t) the 
desired evolution, is a measure of the effect of the unitary 



evolution on encoded quantum information. If pu(t) = 
p(t), then F(t) = 1 (for pure states), and the system 
serves as a perfect quantum memory. The fidelity can be 
expanded as 



t 



for 



Tr[{ Pu (t)p(t)}W] 



(25) 



(26) 



where the superscript (n) denotes n th derivative @, HqI- 
The timescale t^ 1 = for any H5, so we focus on ■ 
This is an estimate of the timescale over which quantum 
error correction (or some other technique , eg . a corrective 
pulse sequence) will have to be applied [271 ]. So, 



= (^|H S |^ 2 -(^|H||^), (27) 



for within a DFS. We have assumed the system begins 
in a pure state in a DFS and that the desired evolution 
satisfies p(t) = p(Q). We are interested in the transfer 
timescale, so we assume that L-£,[pjj(t)] 7^ only for t > 
t 2 . If the states \ip) are eigenstates of Hs, then t^ 1 = 
0. The timescale in Eq. (|27|) also applies to transferring 
quantum information between states within a DFS. 



V. SMALL QUBIT SYSTEMS 

In practical applications and in many theoretical pro- 
posals, the number of qubits that are controlled in order 
to process quantum information is small. Here, we ex- 
amine in detail systems that might be realizable in a lab- 
oratory, and give conditions on fe Q( g that lead to robust 
quantum information storage. We focus on three and four 
qubits because of the possibility that b a p is spatially de- 
pendent (which is the case for dipole-coupled qubits). If 
so, experimental control of b a p can be obtained through 
varying the spatial arrangement of the qubits. 



A. Three qubits 

Three qubits is the smallest number of qubits that sup- 
ports a DF qubit Q . The eigenbasis is given in Fig. [3] 
The qubit is encoded as 



|1>l = 



10), 



lf) : 



^(|010) -|100)), 



V2 



(|011) -|101)), 



1 

i(2|110) 



2|001) + |010) + |100)), 
|101) - |011», 



(28) 
(29) 



where the logical groupings are indicated in the first col- 
umn, the basis used in this section is given in the sec- 
ond column, and the states arc expanded in terms of the 
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10), 



|1111) = |2, 2) 




10 = 1 



|n) = |l,l) |ra) = |l,l) |l)=|l,l) | k) = | 2 .i) 



W = |1,0> |g) = |l,0) |f) = |l,0) |e) = |2,0) 



|d) = |l,-l) |c) = |l,-l) |b) = |l,-l) M = |2,-l) 



W = l|, 



|0000) = |2, -2) 



FIG. 3: DF encoding for three qubits. The states are la- 
belled according to \ J,m.j). The two isolated subspaces are 
circled according to the logical basis. See Eqs. (|28|l and (|29|) 
for details. The splitting of the degeneracy due to the non- 
nearest-neighbor unitary interaction is not included. 



FIG. 4: Eigenbasis for four qubits, labelled \J,mj), with the 
logical DFS explicitly labelled. The splitting of the degen- 
eracy due to the non-nearest-neighbor unitary interaction is 
not included. 



single-particle basis in the third column. Transitions are 
only allowed between states with the same symmetry, i.e., 
the jump operator J does not cause quantum information 
to decay from |e) to |c) or |d), or from |f) to |b) or |d). It 
only acts within the logical groupings. The degeneracy 
for each J is given by Eq. (f2~4"|) . 

The off-diagonal elements of the Hamiltonian 



He 



a, (3=1 



(30) 



where a and (3 label the qubit, in the one-excitation sub- 
space in the collective basis are 

H J S =-^(&23 - 6l3) |d) (C| + ^={b 13 - b 23 ) | C ) (b| 



1 



3V2 



(26 



12 — "13 ~ "23, 



|d)(b| + H.c. (31) 



For dipole-coupled qubits, the coefficient b a p describes 
the dipole-dipole interactions. One possible physical 
system that might satisfy the required conditions are 
nitrogen-vacancy (NV) centres in diamond. These can 
be manufactured in such a manner that the expansion 
detailed in Eqs. Q and (fTU|) adequately describes the 
physical system [281 ) . 

A consequence of [J z , H5] — is that the Hamiltonian 
is block-diagonal in the collective basis. We examine the 
three cases: (i) b a p = 6, (ii) 612 = 623 7^ &13, and (hi) 
b\2 7^ 623 7^= 613. It can be seen immediately for b a p = b 
that H l s = 0, and Eq. (fT5|) is satisfied. It can also be 
seen for 623 = 613 that |c) is not acted upon by H5 and 
is CDF. 

We consider b\ 2 = &23 7^ 613 ■ Introducing the states 
|u) = i |c> - ^ |b) and |v) = & |c) + ± |b), the off- 
diagonal system Hamiltonian becomes 



= ^(633-613)^) 



H.c. 



(32) 



meaning that |u) is not acted upon by H5 and is CDF. 
So, if state c u |u) + c v |v) for |c u | 2 + |c v | 2 = 1 is prepared, 
the second order transfer rate is 

-X ^Ihscl- 1 - (b 13 - 4b 12 )c 2 v } 2 - 2b 2 13 c 2 u 

T2 / 3 

-^(6b 2 12 -4b 12 b 13 + b 2 13 y v , (33) 

which implies that smaller differences between the ele- 
ments of {bap) leads to more robust quantum informa- 
tion storage. The limiting case 623 — > b, 613 — > b recovers 

(i)- 

Consider the most general case: &12 7^ 623 7^ 613 ■ It 
can be seen that there are no stationary states using 



P = Mp, 



(34) 



where M. is the restriction of H5 and Ld to the one- 
excitation subspacc. For the steady state p — 0, and 
there exists a nontrivial solution to Aip = iff det(A^) = 
0. This is calculated to be 

4A 3 

det(M) = (61a - b 23 ) 2 (b 12 - b 13 ) 2 (b 23 - b 13 ) 2 , 

(35) 

where A is the nonzero eigenvalue of (a a p). This shows 
that for det(A^) = 0, two elements of {bap) must be 
equal, which is case (ii) above. Note that there are no 
CDF subspaces for the general case. 



B. Four qubits 

In order to exploit the collective properties of a system 
of qubits, t 2 x gives the timescale over which one would 
expect to be able to encode information without loss. 
However, this timescale might be much faster than the 
timcscales given by the eigenvalues of (a a p). In fact, for 
a generic system of dipole-coupled qubits the timescale 
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for unitary evolution is ~ 10 8 times faster than the decay 
timescale [l6|. This difference in timescales remains for 
a fully rcnormalizcd theory, so for dipolc-couplcd qubits 
DF quantum information will evolve into non-DF states 
much faster than the decay rate of the system. 

This timescale difference implies that one should en- 
code to protect against the effect of H5 before one en- 
codes against nonunitary decohcrence. Here, we give ex- 
plicit examples for different forms of b a p and show how 
one could encode to protect against unitary evolution. 
A significant result presented here is an encoding that 
protects against an arbitrary b a p in the two-excitation 
subspace. The Hamiltonian 



IT 



a,/3=l 



(36) 



causes evolution between states that have the same value 
of mj. So, encoding in the strong-collective DFS spanned 
by {|i> , |j)} where 



|i) =-^(|01)-|10»(|01)-|10)), 



(37) 



|J> 



1 



2 |0011) + 2 1 1100) - |0101) - |1010) 



12 

0110) - |1001>) 



(38) 



will not guarantee stable quantum memory. 

Focusing first on the one-excitation subspace, we notice 
that if &14 = 623 and 613 = &24 then 

H J S = i(6 12 - 634) |a) <d| + (b 2i - b 23 ) |b) <c| + H.c, 

(39) 

for the basis defined in Fig. [4] So, a logical state en- 
coded across {|b) , |c)} is not acted upon by Hg, and 
satisfies Eq. (|15p. Under the same conditions, in the two- 
excitation subspace the off-diagonal terms of the Hamil- 
tonian become 



Hi 



=yf(&23-&24)(|e) (i| + V2|i 
. V2 



j (623 - 612 + b-24 - 634) |e) 01 + H.c, (40) 

which shows that the DFS {|i) , |j)} is coupled to the sym- 
metric state. The states that are not acted upon by Hg 
are not lowest weight states, and so decay through the 
action of Ld • A further condition is required to decouple 
the DF state |i), namely 623 = 624- Then, there are three 
states that satisfy Eq. (fT5]) : |b), |c) and |i), and a logical 
qubit can be encoded using the most convenient states 
for practical applications. Note that we are interested in 
storage times, and do not consider ease of preparation 
and manipulation. See Ref. [l4| for one possible method 
of preparing and manipulating a logical qubit in a col- 
lection of dipole-coupled qubits using globally-addressed 
bichromatic incident fields. 



It should be mentioned that since the coefficient of the 
operator |j) (e| depends on all values of b a p, any pertur- 
bation away from b a p = b will cause information encoded 
in |j) to decohere. This would prevent the use of an ex- 
perimentally controlled (b a p) being used as a single-qubit 
gate on {|i) , |j)}. 

We now relax the constraints on the system Hamilto- 
nian, and consider arbitrary values of b a p. We concen- 
trate on the two-excitation subspace. The Hamiltonian 
can be split into two parts 



H 



f g h 



Jjeij 



(Hg) acts only on {|f) 



(41) 
|h» 



where H^ h 

({|e) , |i) , |j)}). The antisymmetric DF states {|i),|j)} 
are coupled to the symmetric state |e) that undergoes 
nonunitary evolution. The states of interest in the 
single-particle basis are 

|f) = ^(11100) - 10011)), (42) 
Is) = 5(|0H0) + 10101) - 11010) - 11001)), (43) 

l h ) = od 1001 ) _ I 1010 ) + I 0101 ) _ I 0110 ))- ( 44 ) 

In the collective basis, the zero-logical state in the two- 
excitation subspace that, before a jump occurs, is im- 
mune to an arbitrary b a p is 

\0) L = / J_g (n 1 |g) - [h)), (45) 



for fix = -fe(b 



+ b 2i ), n 2 = -±(b 13 + 614 



b 2 3 — ^24)1 and the temporal evolution associated with 
the diagonal terms in H f | h has been absorbed into |g) 
and |h). The one-logical state is then a combination of 
the two remaining eigenstates of H f | h 



1 



( N /fii r Tn£|f) + n 2 |g) + n 1 |h>), 



(n 2 |g) + ni|h) 



xA^+oflf)). 

(46) 

The states {|f) , |g) , |h)} are coupled to each other, but 
are not coupled with states {|e) , |l) , |j)}. Using the en- 
coding given in Eqs. (j45|) and (|46|) means that, up until J 
acts on the two-excitation subspace, the qubit will be im- 
mune to an arbitrary environment induced non-nearest- 
neighbor evolution. For the physical example of dipolc- 
couplcd qubits this is surprising, particularly since it was 
shown in Ref. [29I ] that including the dipole-dipole inter- 
action destroyed any collective-emission behaviour. Also, 
for dipole-coupled qubits in this regime, the timescale 
for unitary evolution is typically ~ 10 s times quicker 
than that for decay, so using the encoding in Eqs. ([43]) 
and (|46p might have immediate benefits to applications. 
We should emphasize that quantum information will de- 
cay within the J = 1 irrcps, so the proposed encoding 
does not support perfect quantum memory for infinite 
time. 
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VI. CONCLUSION 

It is well-known that the detrimental effect of bath- 
induced Hamiltonians is not accounted for by requiring 
L D [p] = 0, for p a DFS Q ■ In this paper ,we stated a con- 
dition, similar to the conditions in Refs. that en- 
sures persistent DF quantum information in the presence 
of a non-nearest-neighbor bath-induced system Hamilto- 
nian. We showed that, in light of an arbitrary system 
Hamiltonian, as the size of the Hilbert space increased, 
the strong-collective DFS is the least suitable subspace 
for quantum information storage. The most suitable 
place to store quantum information in N qubits — if scala- 
bility is important — is probably the subspace with ~ |iV 
excitations. We then gave a timescale over which other 
methods would have to be applied to account for H5. We 



then concentrated on small qubit systems, giving specific 
examples for three and four qubits that we hope will have 
immediate benefit to applications. A particularly inter- 
esting result for four qubits was the encoding that elimi- 
nates the need to correct for H5 while the qubit remains 
in the two-excitation subspace. 
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